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I. INTRODUCTION 



The photoemission and magnetic properties of cuprate perovskites have been exten- 
sively studied during the last few years, both because of their unusual behaviour and in 
the hope that they might provide insight into the physical origin of high-temperature 
superconductivity. Among these properties the pseudogap observed in photoemission 
and the magnetic pseudogap revealed in the static susceptibility and in the spin- 
lattice relaxation rate of normal-state underdoped cuprates have attracted consid- 
erable attention. A number of different approaches were suggested for the description 
of the pseudogaps. In particular, in view of the similarity in symmetry and size of 
the photoemission pseudogap with the superconducting gap, in Refs. this pseu- 

dogap was connected with the superconducting fluctuations existing above T c . This 
idea was based on earlier theoretical results of Refs. |],|TIJ]. Another point of view 



was suggested in Ref. [11| where the energy spectrum of the two-dimensional (2D) t-J 
model was shown to have a peculiarity which is similar by its properties to the pho- 
toemission pseudogap. In accord with Ref. the pseudogap is a consequence of a 
specific dispersion of the strongly correlated electron system at moderate doping and 
is not connected with superconducting fluctuations. The discussion of the magnetic 
pseudogap in doped cuprates is mainly based on scaling arguments |l2Hl4f and on the 
idea of real-space pairing |T5f . 

In this paper we describe the energy spectrum, including the photoemission pseu- 
dogap, and the magnetic properties demonstrating the magnetic pseudogap in a unified 
approach based on the 2D t-J model widely used for the description of Cu0 2 planes 
of cuprate perovskites (the extensive literature on this model is reviewed in Ref. (jHJ ) . 
For the consideration of the paramagnetic state we extend the spin-wave theory with 



1 



the constraint of zero staggered magnetization, developed for the Heisenberg model 



17|,18|, to the t-J model. The spectrum is determined by solving numerically a set of 



self-energy equations for hole and magnon Green's functions in the self-consistent Born 



approximation with account of the constraint ||19|| . The constraint can be fulfilled in 
the ranges of hole concentrations 0.02 < x ^ 0.17 and temperatures T < 100 K. In this 
region the obtained hole spectrum differs from a conventional metallic spectrum as the 
quasiparticle weights of states are less than 1 and change with the hole concentration. 
This leads to the violation of Luttinger's theorem. 

The hole spectrum consists of two essentially different parts: a persistent portion 
of the narrow spin-polaron band, which is typical for the low-concentration (x < 0.02) 
spectrum, and a wider part appearing from x ~ 0.04 which is characterized by the 
energy parameter t, the hopping constant. The former part provides the most intensive 
features in the hole spectral function near the Fermi level. For x < 0.17 this part, 
which is pinned to the Fermi level near (±7r/2, ±7r/2), bends upwards (in the hole 
picture) on approaching (±7r,0), (0, ±7r). In the hole spectral function the crossing 
of the Fermi level by the second, wider part is completely lost within the foot of a 
more intensive spin-polaron peak in this region of the Brillouin zone. This looks like 
the disappearance of a part of the Fermi surface and the opening of a pseudogap 
near (±7r, 0), (0, ±7r). Obtained size, symmetry and concentration dependence of the 
pseudogap are in agreement with photoemission data in Bi 2 Sr 2 CaCu 2 08 +( 5 (Bi2212) 
11- 

Another peculiarity of the calculated hole spectrum is an extended saddle point 
near (0, ir). By the energy position and by the extension in the Brillouin zone our 
results reproduce well the analogous feature of the photoemission spectra. In the 
case of optimal doping, x ~ 0.17, this peculiarity leads to a 2D Fermi surface in the 
considered 2D system. 

A gap in the magnon spectrum of the undoped antiferromagnet appears as a conse- 
quence of the constraint of zero staggered magnetization ||i~7| , |18fl . Starting from x ~ 0.02 
this gap is filled by overdamped magnons. Their increased damping is the consequence 
of the hole-magnon interaction and at T = it indicates the destruction of the long- 



range antiferromagnetic order by holes p0| . The arising pseudogap leads to the de- 



crease of the static spin susceptibility and the spin-lattice relaxation rate (TiT) -1 with 
decreasing temperature, as observed experimentally. This behavior is typical for the 
quantum disordered regime. Calculated values of these quantities are in qualitative and 
in some cases in quantitative agreement with experiment in underdoped YBa 2 Cu30 6+< 5. 
In some articles the photoemission and magnetic pseudogaps are identified (see, e.g., 
Ref. ||). In our opinion these pseudogaps are two different, unconnected peculiarities 
of the spectra of two well-distinguishable subsystems — holes and magnons. 

The outline of the paper is as follows. In Sec. II we discuss the derivation of the 
t-J Hamiltonian from the extended Hubbard Hamiltonian which is supposed to give 
a realistic description of Cu0 2 planes. This allows us to rewrite spin operators of 
oxygen and copper ions in terms of operators of the t-J model. The extension of the 
modified spin- wave approximation with zero staggered magnetization on the t-J model 
is discussed in Sec. III. The self-energy equations for the hole and magnon Green's 
functions and the hole and magnon contributions to the magnetic susceptibility in 
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terms of these functions are written out in Sec. IV. Numerical results on the spectrum 
are considered in Sec. V. The pseudogap in the hole spectrum and magnetic properties 
are discussed in Sec. VI and VII, respectively. Finally our conclusions are given in 
Sec. VIII. 



II. THE EFFECTIVE HAMILTONIAN 



The extended Hubbard model [21] is widely used for the description of Cu0 2 planes 



of cuprate perovskites. The Hamiltonian of the model can be written in the form 

H = Y,H m + 2t pd X a J2 (4u70m+a,«r + H.c.), 

m macr . ^ . 

H m = C/n m;+ in m _i + A^0 ma m(T + 2t pd X ^(d^^m* + H.c), 

cr cr 

where c? mo . is the creation operator of electrons in the 3d x 2_ y 2 orbitals of copper at 
the plane site m with the spin a — ±1, (jy ma . is the Fourier transform of the operator 
^ko- = (Pk/2y/N) Ema ex p( — ^ m )P m +a/2 cr constructed from the creation operators of 
electrons in the 2p a orbitals of oxygen p m+a / 2 a - Complementary linear combinations of 
these operators, which do not hybridize with the 3d x 2_ y 2 copper orbitals, are omitted 
in Eq. ([]]) because their energy is much higher. In Eq. ([[]), a = (±a, 0), (0, ±a) where a 
is the in-plane copper distance which is taken as the unit of length, /3k = {1 + [cos(k x ) + 
cos(k y )]/2}~ 1 / 2 , N is the number of sites; U, A, and t p d are the Hubbard repulsion on 
copper, the Cu-0 promotion energy and hybridization, respectively, n mu = rf mCT d m(7 , 
A m = N' 1 Xikexp^km)^ 1 , A ~ 0.96, A a ~ 0.14. Other components of A m are small 
and the respective terms are omitted in Eq. (|l]). 

The splitting of the Hamiltonian into the one- and two-site parts in Eq. (|I|) pro- 
vides a good starting point for the perturbation theory, because for parameters [23| of 



La 2 Cu0 4 these two parts are characterized by energies differing by one order of mag- 
nitude. Notice that the frequently used perturbation expansion in powers of t p d does 
not work for these parameters, as the hybridization is actually not small in comparison 
with other energies f22| . |2~4}| . The pref actor A a ~ 0.1 in Eq. (jl|) allows one to overcome 



this difficulty The zero-order, one-site part ^ m H m of the Hamiltonian has two sets 
of states corresponding to unoccupied and occupied site states of the t-J model 



l m ) — "y|(0m,+l<4n-l — 0m,-l^m,+l) + C 220m,-10m,+l + C 23^L,-1^L,+1 
|ma) = (c310L,-10L,+l^L + C320L rf L,-l^L,+l)l^m), 



\v m ), 

(2) 



where \v m ) is the site vacuum state and the coefficients are obtained in the course of 
the diagonalization of H m . For a given number of holes crystal states \q) constructed as 
products of site states (@) form the degenerate ground states of the zero-order Hamil- 
tonian H = J2 m Hm- These ground states are separated by a finite gap of the order 
of min(A, t p d) from excited states. In such conditions one can use the operator form of 
the perturbation theory G5| to obtain an effective Hamiltonian acting in the subspace 
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of the low- lying states \q). Up to the terms of the second order in the perturbation Hi 
[the two-site part of Hamiltonian (Q)] this effective Hamiltonian reads 

H eS = V\H X - - V)(H - E )-\l - V)Hi]V, 

where V = J2 q \q)(q\ an d Eo is the eigenenergy of H for the ground states \q). Using 
Eqs. (JXp and (|2]) we obtain the effective t-J Hamiltonian [^ 

H e fi = tY^ |m + a,<7}(m + a||m)(mo-| + ^S m S m+a , (3) 

macr * ma 

where S* m = S^ + icrS^ = |mcr)(m, — a\, 5* m = X^(^/2)|mcx)(m<7| are the components 
of the spin operator S m , t and J are the effective hopping and superexchange constants 
which are expressed in terms of the parameters of the extended Hubbard Hamiltonian 
[see Ref. p2]| ; in Eq. @, we omitted three-site (t') terms which give only small correc- 
tions in the considered case J/t <C 1] . Using parameters of Ref. |2|j we estimated the 



ratio J/t to lie in the range 0.1-0.5. In the following discussion we use J/t = 0.2 and 
t = 0.5 eV. 

The Zeeman term of the Hamiltonian can be written in the form 

Hz = 2fi B ^2 S mH m + fiB ^2 S m+a/2H m+a / 2 , (4) 
m ma 

where fiB is the Bohr magneton, H m is the applied magnetic field, and s m and s m+a / 2 
are composed of 

s' 7 = r/t d h z = V —d^ d 

cr Z 

cr _ f Z _ V"^ ® t 

S m+a/2 — Pm+a/2,crPm+a/2,-cr, S m+a/2 _ ^Pm+a/2,aPm+a./2,a , 

cr Z 

respectively. Using notations of Eq. ([!]) the oxygen spin operators can be approximately 
rewritten as 



/3 2 

S m+a/2 = -rifimo + 4+a,a) (0m,-cr + 0m+a,-cr 



Sm+a/2 = T Yl + 0m+a,cr)(0mcr + 0m+a,<x), 

where /3 = /? + /3 a ~ 0.96, j3 m being the Fourier transform of (3^. In the basis of states 
(0) the spin operators read 

Sm+a/2 = C2iC3 1 /3 2 ^^(|ma)(m||m + a)(m + a,or| 

+ |m + a, cr) (m + a||m)(mcr|), (6) 

Sin = C3lEljl m(T )( m(J l> 
cr Z 

where in the oxygen operator terms containing small coefficients C32, C23 were omitted. 
In the following discussion we set C21 = C31 = j3 — 1. 
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III. THE SPIN- WAVE APPROXIMATION 

FOR ZERO STAGGERED MAGNETIZATION 



Hamiltonian can be essentially simplified with the use of the spin-wave approx- 



imation. As known |26|, at low temperatures and hole concentrations x < 0.02 the 
CuC>2 planes are antiferromagnetically ordered. For larger x this long-range ordering 
is destroyed. To describe low-lying magnetic excitations and their interaction with 



holes in this case we use the version of the spin- wave theory formulated in Refs. fT7|Ti"8 l 



for the Heisenberg antiferromagnet with zero staggered magnetization. As shown in 
Refs. [P~7| . p~S[ ] , in the absence of holes this approach reproduces results obtained in 



Refs. [p7| , p8H with the mean-field Schwinger boson and renormalization group theories 
and it is remarkably accurate, as follows from the comparison with exact diagonal- 
ization and Monte Carlo results. Here we reformulate this approach to simplify the 
inclusion of holes in the Hamiltonian. 



We use the Holstein-Primakoff transformation |25[ to introduce boson operators of 
spin waves 6 m , 

SL = e mm Q - &L^m) , S° m = P^ m b m + P-'bl^, (7) 

where II = (tt, 7r), P m = [1 + aexp(iIIm)]/2, and <^ m = (1 - & m & m) 1/2 - In Eq. (0), 
the factors exp(illm) and P^ 7 are introduced to account for alternating directions 
of spins in the classical Neel state which is used as the reference state in the spin- 
wave approximation. On substituting Eq. (|7|) into the Heisenberg part of Hamiltonian 
(H), expanding ip m and keeping terms up to the quartic order, we use the mean- field 
approximation in these latter terms 



J ^ „ „ JN 



H H = 2 ^ S m+aS m « — - J (b b a ) 4 b\j) m + - (^m+a^L + b m+a b 



1 



2 



•t , „ 

/ m+a L 'm 



(8) 



where angular brackets denote averaging over the grand canonical ensemble and the 
four correlations {b b a ) are supposed to be equal. On deriving Eq. @ we took into 
account the condition 

(blb m ) = \ (9a) 
which follows from the constraint of zero staggered magnetization, 

E^X = or £&U»=y (9b) 

and ensures zero site magnetization, (S^) = 0. To account for this constraint in 
the subsequent consideration we add the term 2Jv(bab a ) Z^m^L^m with the Lagrange 
multiplier v to Hamiltonian @. 

The deviation of the quartic terms from their mean-field value in Eq. (|8|) describes 
the magnon-magnon interaction which leads to the magnon damping |29]| . In the 
t-J model there is another mechanism of the magnon damping connected with the 
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hole-magnon interaction. Estimations p0| , p9| show that in the considered range of 
hole concentrations this latter interaction gives the main contribution to the damping. 
Therefore we do not consider the magnon-magnon interaction below. 

The resulting Hamiltonian is diagonalized by the unitary transformation 



U 



exp 



Z k<7 



&L&U,- 



(10) 



with «k = l n [(l+^7k)/ (1 — 7 77k)]/4, r\ = 2/ (2 — v), 7k = X)a ex P(*ka)/4, and the primed 
sum sign indicates that the summation is restricted to the magnetic Brillouin zone 
which is half as large as the usual one. In Eq. fllPp, 6 k(T : 



/2/iVE m exp(-ikm)6 m P- 
where due to the projector the summation is performed over one sublattice. As a 
result, we obtain 



4J 



7k 



{b-k-ab^u - r]7 k (bipbka) 



'ID 



where we omitted unessential constant terms, and in (bob a ) for the following discussion 
we keep the anomalous correlation (6-k,-o-^ko-)c which is nonzero at i ^ 0. The 
subscript U means that the averaging is performed with the Hamiltonian transformed 
with operator fllPf) . In the absence of holes we have (6 ko .6ko-)t/ = [exp(a; k /T) — 1] _1 
and (6-k,-cr&ko-)!7 = where T is the temperature in energy units, while for x > the 
correlations are calculated from the magnon Green's function. Condition (|9]a) which 
determines rj in Eq. ( [TT1) acquires the form 



N 



1 



=— r ((&L & ka)c/ + ^ - ?77k< 
- ^ 2 7k V A 



^-k-abka)u 



1. 



(12) 



Analogous equations (without (&-k,-cr&k(r)i/) for the magnon spectrum of the Heisen- 
berg antiferromagnet without holes were obtained in a somewhat different manner in 
Refs. ]T7| , [r8| ]. As shown in these works, in states without long-range antiferromagnetic 
ordering and for finite lattices one obtains rj < 1 which introduces a gap in the magnon 
spectrum (|TT|) near the points (0, 0) and (ir, n) of the Brillouin zone. 

The reference state of the spin-wave approximation discussed above is the classical 
Neel state | J\f) . Other states are described via the reference state and magnon creation 
operators determined for this state. The introduction of holes in this picture leads to 
two possibilities for the hole movement: there is a magnon or there is no magnon on 
a site which a hole jumps to. Both these possibilities are described by the following 
term of the Hamiltonian: 



H t = tJ2 h m hln +a (b m+a + 6 m ), 



(13) 



which corresponds to the first term in Eq. (|3|). In Eq. (|TB|), = J2a -P m l m ) ( mo "| is 
the hole creation operator in the Neel state \J\f) = Y[ m iJ2a -P m l mcr ))- 
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In spite of the clear physical meaning of the constraint of zero staggered magne- 
tization, let us discuss the approximations made from a somewhat different point of 
view. Notice that in the case of short-range antiferromagnetic order characterized by 
the spin correlation length £ ^> a one can use the usual spin-wave approximation in 
any crystal region with a linear size of the order of £. A Neel state with some local 
spin quantization axis is used as the reference state for this spin-wave approximation. 
Additionally one should take into account the finite-size effect due to a finite value of £. 
This is done by applying condition (0). The value of rj determined by this condition is 
directly connected with £ JHJ (recall that r\ < 1 both for a finite £ and a finite lattice). 
Notice also that only sites neighboring to the hole site are involved in the processes 
described by Eq. fll3|) . The hole-magnon interaction is of the short-range type. This is 
the reason why in the case of short-range antiferromagnetic order with £ 3> a the hole 
hopping term (|i~3"D looks exactly like the same term for the long-range order |]30| , |3T 



There is no "bare-hopping" (without magnon operators) term in Eq. ( |T3"| ) because in 
the case £ ^> a the direction of the spin quantization axis is practically the same on 
sites involved in a hole jump. 

After unitary transformation (llOj) the total Hamiltonian reads 



H = U^H t U + H H = E W4A-k',-A' CT + H.c.) + £^k&kA CT , (14) 



kk'cr kcr 



where = — 4ty'2/iV(7 k _ k /Mk' + 7k^k') and = cosh(ak), t>k = — sinh(«k)- On 
carrying out the unitary transformation we neglected the noncommutativity of hole 
and magnon operators which is justified at least for x ^5 0.1 by comparison with exact 
diagonalization results |32|]. At r] — 1 Eq. fll4|) reduces to the Hamiltonian obtained 



in Refs. fl30| , |3T|j for the case of the long-range antiferromagnetic order. The number of 
magnons is not conserved by Hamiltonian ([14]) and at x > the anomalous correlation 
{b-k-abk^u in Eqs. (|i~T|) and (|12|) is nonzero. 

Notice that the used spin-wave approximation is not rotationally invariant and the 
correlations {S^S^} are zero 1 17|,18| . Therefore only the z components of spin operators 



are considered in the following discussion. 



IV. SUSCEPTIBILITY AND SELF-ENERGY EQUATIONS 



In the new notations the spin components (||) acquire the form 



Ilm 



- - b ] b 

2 u m u m 



i i 'm+a/2 



illm 



32 



h h ] 



i n ' ' m \ a I ^m+a 



bl 



H.c. 



• (15) 



and s* , /2 give contributions to the magnetization 



'm+a/2 



-iq(m+a'/2) 



from magnons and holes, respectively. Here a' = (a, 0), (0, a). The susceptibility is 
determined by the equation 



M z . .\l 

qr> — q 
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where M^ T = exp(i#r)M* exp(-iftr), H = H t + H H - /jM, M = Eka4<A<x, and 
is the hole chemical potential. The susceptibility can be calculated using the sim- 
plest decoupling, as the subsequent terms of the perturbation series for the respective 
Matsubara Green's function are proportional to powers of the small hole concentration 
and are further decreased by rapidly oscillating coefficients. In this approximation the 
magnon and hole contributions (indicated by subscripts m and h, respectively) in the 
susceptibility can be written as 

Imx^(q^) = ^Ef — [n B {v)-n B {v + u)\ 

k 00 71 

x [^(k, z/)ATi(k + k, v + u) + K 2 (k, v)K 2 (k + n,u + cu)}, 
lm X z h (quj) = II —^-[ n F(. u ) ~ u f( u + u)][n B (v') - n B (i/ - u)] 



kk' _uu 7r 



x Im G(ku) Im G(k + k' - q, v + 1/ - u) (16) 



x 



Im D n (kV) 7 2 (k - |) - Im D u (k', -v') 7 2 (k + k' - |) 



-2]mr; J2 (k / //)-(k-|J 7 (k + k , -| 



re 

Rex(qc^) =V / 



Im x(q a; ) 



'-00 7T V — W 

where n B = [exp(a;/T) — 7ip = [exp(u/T) + n = q — II, 

Ki(kco) = w k Im Dn(kuj) + 2-Uk^kIni Di 2 (kuj) — v^Im -Dn(k, —a;), 

if 2 (ku;) = u k f k [ImDn(ka;) - ImDn(k, —u)] + (u 2 , + Vk)ImD 12 (ka;), 

and Dij(\tLd), G(ku) are the Fourier transforms of the magnon and hole Green's func- 
tions 

D n (kt) = -z0(*K[M*),&L])tf, ^(kt) = -i0(*)([&- k ,- CT (i)A<x]>t/, 

G(kt) = -^(t)([/i kCT (t),4 (T ]) l/ , 

^ko-(^) = exp[i(7Y — /xA/")t]6k<T ex P[ — i(H — H>N)t\. These Green's functions do not depend 
on the spin index. In Eq. (|T6|) , summations over wave vectors are performed over the 
full Brillouin zone. In the second magnetic zone, which together with the first magnetic 
zone forms this full zone, Du(ku) and G(ku) repeat periodically their values in the 
first zone, while D 12 (kuj) changes sign. 

These Green's functions are determined from the set of self-energy equations. The 
hole S and magnon fin, II12 self-energies are described by the following diagrams: 



= O ( • + O 4 • + C ( • + G <J • . 



where solid and dashed lines correspond to hole and magnon Green's functions, open 
and filled circles are bare and full vertices. Dashed lines with two oppositely directed 
arrows correspond to the anomalous magnon Green's functions -D 12 (kt) and _D 2 i(kt) = 
— i@{t)([bl_ a (t), fe! k _ a \)u- The first correction to the bare vertex g kk / is exactly zero due 



to the impossibility to coordinate spin indices in this correction |33|. This suggests the 
use of the Born approximation in which the full vertices are substituted with the bare 
ones. In this approximation the real-frequency self-energy equations read 



G{ku) = [u + n - ^(Mn 1 
Im£(ko;) - 



2 E / - 4lm5 U (kV) - sL k , _ k JmD n (-k', 

j^, J — oc 71 L 



-0J 



+25'kk'5 , k-k',-k'ImDi 2 (k / u; / ) [n B {u') + n F (uo' - w)]ImG(k — k',u — 10') 
duo' ImS(W) 



oo 71 UJ' — UJ 

R*(k,-u) 



Re£(ku;) = V 
Dn(kw) 

U{KOJ)n^{K, - 

i roc 
Imn n (ka;) = 2^^ k 



R(ku)R*(k,-u) -nf 2 (k,a; 

-oo 71 

00 du' 



D l2 {ku) 



ni 2 (k,a;) 



(17) 



R(ku)R*(k,-u)-Ill 2 (k,u) 
Im G(k'u') Im G(k' - k, u' - u) \n F {uj') - n F (u' - u) 



Imn 12 (ka;) =2^ <? k ^_ k ,_ k / — ImG(kV) ImG(k' - k, J 

J— oo 71 



UJ 



x [n F (u') — n F (u' — uj)], 

00 du' Imlly(W) 



Rell^(kw) = 7> 



00 7T 



UJ — UJ 



where R(ku) = uj — u; k — IIii(ka;). Equations fll7D and ( |12D form a self-consistent set 
and can be solved iteratively. The calculation procedure is the following: for given 
values of /1 and T some starting value rj < 1 was selected; after the convergence of 



300 



200 



T, K 




Fig. 1. The curve encloses the region where 
condition (fl^) can be fulfilled. 
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the iterations condition (|12|) is checked and rj is appropriately changed for the next 
iteration cycle, until the condition is fulfilled with the accuracy of 10~ 3 (some other 
details of the calculation procedure can be found in Ref. [0). In the calculations 
a 20x20 lattice was used. Green's functions were computed on a mesh of frequency 
points equally spaced with the step Au; « 0.022t « 11 meV. We found that condition 



(|12|) can be satisfied only in a limited region of the T-x plane, namely below the curve 



shown in Fig. 1. 



V. THE ENERGY SPECTRUM 



The evolution of the calculated hole spectral function v4(ku) = — Im G(ku) with 
the concentration 



2 poo 

x — / duj np(uj)lmG(kuj) 

nN , J-oo 



k 

is shown in Fig. 2 for two points of the Brillouin zone. We use the hole picture where 




states below the Fermi level u = are filled by holes. Since our description is based on 
the spin-wave approximation with the magnetic Brillouin zone which is twice smaller 
than the full Brillouin zone, the spectral functions in the points k and k + (jr, n) are 
identical. For long-range antiferromagnetic ordering when these points are equivalent 
such description does not lead to any loss of information. In the case of short-range 
order the use of the smaller Brillouin zone leads to a somewhat coarsened description. 
The points k and k+ (it, tt) are no longer equivalent — the respective spectral functions 
contain similar maxima which however have essentially different intensities and widths 
P5fl . In our approach, where the maxima in the two points appear together in one 
spectral function for the momentum in the first magnetic Brillouin zone, we cannot 
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determine whether a selected maximum is more intensive in the first or in the second 
magnetic Brillouin zone (which form the full zone). In other words we cannot determine 
where is the usual band and where is the shadow band . However, we can distribute 
the maxima between the two magnetic zones by drawing additional information from 
experiment as we shall do in the following discussion. 

At x < 0.04 the spectra contain series of peaks and the hole-magnon scattering 
continuum (see curves 1 in Fig. 2). With increasing x only the lowest and most in- 
tensive of these peaks is retained in the spectrum for wave vectors near the boundary 
of the magnetic Brillouin zone [everywhere in this region the spectrum is similar to 
that for k = (0, 7r)]. This peak corresponds to the so-called spin-polaron band in the 
hole spectrum. Other peaks are washed away, forming a broad dispersive maximum 
above the Fermi level for momenta in the central part of the magnetic Brillouin zone. 
Simultaneously a maximum develops below the Fermi level (compare curves 2 and 3 
in Fig. 2). As seen in Fig. 3, this maximum possesses dispersion and is comparatively 
narrow near the Fermi level. 

The energy positions of the spectral maxima discussed above as functions of mo- 
mentum are shown in Fig. 4 for heavily underdoped and moderately doped cases [as 
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Fig. 3. The hole spectral function for k = 
(0.57r,7r) (a), (0.4vr,7r) (b), and (0.3vr,7r) (c). 
x = 0.133, T = 0. 
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Fig. 4. Energy vs. momentum relationships calculated for x = 0.021, T = (a) and x = 0.121, 
T = 0.02i = 116 K (b, solid line). In part (b), circles indicate positions of quasiparticle peaks 
in the photoemission experiment [4C] carried out at T = 100 K in a Bi2212 crystal with 
T r = 85 K. 



discussed above, we have used the experimental fact that hole pockets are positioned 
around (±7r, ±7r) to relate the portion of the energy band below the Fermi level to this 
part of the full Brillouin zone]. Figure 4a demonstrates the well-known spin-polaron 
band discussed in a large number of works devoted to the case of low doping (see, 
e.g., [j^ , [3^ ,[36H39| ) . The width of this band is of the order of J, i.e. much less than 
the hole hopping constant t. This is connected with the above discussed fact that 
in the rigorous antiferromagnetic order the hole movement requires the emission and 
absorption of magnons. In such conditions the slower subsystem — magnons — will 
determine the bandwidth of the combined quasiparticle. As mentioned, a part of the 
spin-polaron band is retained near the boundary of the magnetic Brillouin zone for 
moderate doping. As seen in Fig. 4b, this part retains, with some distortion, main 
features of the spin-polaron band in this region. Among these features is the large 
nearly flat region which is visible around (0, 7r) in Fig. 4b. In this figure we have 
compared the calculated energy band with the normal-state photoemission data 
in Bi2212 with T c = 85 K which corresponds to the hole concentration in the range 
0.11 — 0.15. In Ref. |4(J the energy position of the flat region — the so-called extended 
van Hove singularity — is estimated to be ±30 — 50 meV relative to the Fermi level. We 
positioned the experimental flat region at +20 meV (in the hole picture) , in accordance 
with later more exact measurements |J. As seen in Fig. 4b, with respect to both the 
energy position and the extension in the Brillouin zone the calculations reproduce well 
the experimental van Hove singularity. Notice that also the general shape of the band, 
as seen in photoemission (which tests the region u > 0), is well reproduced by the 
calculations. That the flat portion of the spin-polaron band can correspond to the 
photoemission extended saddle point was apparently first indicated in Ref. JIT] . 

As follows from Fig. 4, considerable changes occur in the hole band shape on moving 
from light to moderate dopings, as was first indicated in Refs. [B^,B4l] • The narrow, with 
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the width of the order of J, spin-polaron band is transformed to a much wider band 
characterized by the energy parameter t. The general shape of this band resembles 
a considerably distorted 2D nearest-neighbor band produced by the kinetic term of 
the t-J Hamiltonian (^). Analogous changes in the band shape are observed in the 
photoemission of Bi2212 0. These changes point to a certain weakening of correlations 
for x > 0.04, however, some features of the strongly correlated spectrum are retained: 
as mentioned, near the Fermi level around the boundary of the magnetic Brillouin zone 
the spectrum contains a persistent part of the spin-polaron band and widths of spectral 
maxima grow steeply with distance from the Fermi level. 

These changes in the hole spectrum are connected with changes in the magnon sub- 
system. For T = starting from x ~ 0.02 in the central part of the magnetic Brillouin 
zone magnons become overdamped. The overdamped magnons manifest themselves 
in a perceptible intensity of the magnon spectral function Bu(kuj) = — Im Dn(kuj) in 
the nearest vicinity and on both sides of u = (see Fig. 5a where together with the 
structure corresponding to an overdamped magnon the maximum of a usual magnon 
is also visible; notice that overdamped magnons do not form maxima in Bu because 
this quantity changes sign at the central frequency u = of these magnons). The 
appearance of the overdamped magnons points to the destruction of the long-range 
antiferromagnetic order by holes [20,|2| which, in contrast to the destruction due to 
thermal fluctuations, occurs also at zero temperature. Due to a finite intensity in 
Bu, produced by these magnons in the range u < 0, the magnon occupation num- 
ber nk = — J^ 00 duJ7r~ 1 nB{uj)Bii(kuj) is finite at T = which leads to a finite zero- 
temperature correlation length £ in the spin correlation function (sfs^) (in the con- 
sidered finite lattice this phase transition is smeared; fortunately for relevant small x 
an analytic consideration for an infinite lattice is possible The overdamped 

magnons can be identified with relaxational modes describing relative rotations of mag- 
netic quantization axes in regions of size £. In the hole spectrum the destruction of 
the long-range antiferromagnetic order manifests itself in the mentioned change of the 
shape and characteristic energy of the spectrum from J to t at x ~ 0.04. Qualitatively 
this can be understood in the following way: after the destruction of the long-range 
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Fig. 5. The magnon spectral function Bu(kuj) for k = (0, 7r/5), T = 0, x = 0.059 (a), 
k = (0, 0), T = 58 K, x = 0.1 (b), and k = (vr/2, vr/2), T = 116 K, x = 0.172 (c). 
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order holes can move without introducing additional disorder in the magnon subsystem 
and, as a result, the larger characteristic energy, t, reveals itself in the spectrum. 

Apparently nearly simultaneously with the appearance of overdamped magnons 77 
becomes finite and a gap opens in the magnon spectrum. The low-frequency over- 
damped magnons transform this gap to a pseudogap. This pseudogap is most evident 
in -Bn(k = 0, uj) by the asymmetry of the spectrum around u — (Fig. 5b). In spite 
of the appearance of overdamped magnons, even at x ~ 0.17 and T « 100 K usual 
magnons with essentially softened frequencies and increased damping are retained at 
the periphery of the magnetic Brillouin zone indicating persistent antiferromagnetic 
correlations (see Fig. 5c). 

Let us return to the hole spectrum. The hole Fermi surfaces calculated for moderate 
dopings are shown in Figs. 6a and 6b. In part (a), line segments along the boundary of 
the magnetic Brillouin zone are connected with the spin-polaron band. After touching 
the Fermi level (which occurs at x ~ 0.01) the bottom of this band remains pinned to 
the Fermi level and the band flattens with increasing x. Originating from work |43| 



it is widely believed that for small hole concentrations the Fermi surface consists of 
small hole pockets around (±7r/2, ±7r/2). Our calculations do not support this point of 
view. The Fermi level does not cross the bottom of the spin-polaron band but is rather 
pinned to this bottom and the mentioned hole pockets do not arise [notice however 
that the line segments near (±7r/2, ±vt/2) in Fig. 6 can have a finite width provided 
that it is much less than our momentum resolution 7r/10]. 




Fig. 6. The calculated hole Fermi surface for T « 100 K, 0.07 < x < 0.17 (a), x = 0.172 (b), 
and experimental Fermi surface for several cuprates (c, from Ref. fi4|| ). Dashed lines in (a) 
and hatched regions in (b) indicate hidden and two-dimensional parts of the surface. The 
points X, Y, and M correspond to k = (n,0), (0, — vr), and (vr,7r), respectively. 



In Fig. 6a, dashed curves indicate the crossings of the Fermi level by the part of 
the energy band which arises with growing x below the bottom of the spin-polaron 
band (see Figs. 2-4). This part forms hole pockets around (±7r, ±7r) and determines 
the hole concentration (fHf ). The pockets start to form from small x and their extent 
in the Brillouin zone grows with x; however up to x ~ 0.07 the pockets do not reach 
the Fermi level — all their states are positioned below the Fermi level. Thus, for 
x < 0.07 the Fermi surface consists of only the line segments along the boundary 
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of the magnetic Brillouin zone and no closed Fermi surface exists. For larger hole 
concentrations the hole pockets reach the Fermi level and the crossings depicted by 
the dashed curves in Fig. 6a arise. However, as seen in Fig. 3, just at this crossing 
the weaker maximum corresponding to the pocket is completely lost within the foot of 
a much more intensive spin-polaron peak which is located somewhat above the Fermi 
level for the temperature, hole concentration and wave vectors of Fig. 3. Thus, the 
dashed curves in Fig. 6a correspond to hidden parts of the Fermi surface. If like in 
photoemission experiments the spectrum is not tested below the Fermi level (in the hole 
picture), the respective crossings reveal themselves only as a finite spectral intensity 
at the Fermi level. 

As mentioned, with growing x the persistent part of the spin-polaron band flattens. 
As a consequence, the extended saddle points approach the Fermi level and at x = 0.172 
they fall on it. The Fermi surface of the considered two-dimensional fermions becomes 
two-dimensional, as shown in Fig. 6b. The mentioned hole concentration is close to 
optimal doping which corresponds to the highest T c . As observed in photoemission of 
Bi2212 0, in this case the saddle points do lie, within the experimental accuracy, on 
the Fermi level. Fig. 6c demonstrates Fermi surfaces deduced from photoemission in 
Nd 2 _ x Ce !E Cu04 (NCCO), Bi2212, and Bi 2 (Sro.9 7 Pr .o3) 2 Cu0 6+( 5 (Bi2201). Contrasting 
this figure with Figs. 6a and 6b, we conclude that our calculations reproduce satisfac- 
torily the main features of the experimental Fermi surface in crystals with significantly 
different hole concentrations. 

For Fig. 6a we indicated the concentration range where the Fermi surface remains 
practically unchanged. In particular, it means that the size of the (±7r, ±7r) hole 
pockets varies only slightly in the concentration range 0.07 < x < 0.17. This result 
of our calculations agrees with experiment [|TJ and points to the substantial violation 



of Luttinger's theorem ||45j| . In accord with this theorem the area enclosed by the 
Fermi surface in the Brillouin zone should vary linearly with x. As seen in Fig. 2, 
quasiparticle weights of states (spectral intensities) which form the hole pockets grow 
with increasing x. It is this growth of quasiparticle weights, rather than the growing 
size of these pockets, which leads to the increase of the hole concentration in Eq. (p!8|). 

As follows from the above discussion, the energy spectrum of the considered 
fermions differs essentially from the Fermi liquid behavior of conventional metals. To 
this we can add the linear, rather than quadratic, frequency dependences of the hole 
decay widths near the Fermi level in the case of moderate doping []3~5 |. This result 



resembles the marginal Fermi liquid [f46| |, however, in contrast to this concept our cal- 
culated Im£(ci;) has markedly different slopes below and above the Fermi level and a 
strong momentum dependence. 



VI. PSEUDOGAP IN THE HOLE SPECTRUM 



As follows from Fig. 3, the spin-polaron peak is the most prominent feature of the 
hole spectrum near the Fermi level for wave vectors in the vicinity of the Fermi surface. 
The hole spectral function is directly related to the photoemission spectrum. Hence the 
spin-polaron peak will determine the position of the leading edge of the photoemission 
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Fig. 7. The hole spectral function for wave vec- 
tors on the Fermi surface. Curves from top to 
bottom correspond to k = (0.27T, 7r), (0.2-7T, 0.97r), 
(0.27r,0.8vr), (0.3vr,0.7vr), (0.4tt, 0.6vr), and (0.5vr, 
0.5vr), respectively. T = 116 K, x = 0.121. 
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spectrum. For several wave vectors on the Fermi surface the hole spectral function is 
shown in Fig. 7 for the underdoped case. As seen from this figure, the spin-polaron 
maximum, which lies on the Fermi level near (7r/2,7r/2), is shifted upwards on ap- 
proaching (it/ 5, it). Recall that the part of the Fermi surface near (7r/5, tc) is connected 
with another, broader and weaker maximum which is completely lost at the foot of 
the spin-polaron peak on crossing the Fermi level (see Fig. 3). As a consequence, the 
situation in Fig. 7 looks like a part of the Fermi surface disappears and a gap opens 
between the hole energy band and the Fermi level near (7r/5,7r). Due to the men- 
tioned hidden crossing of the Fermi level in this point there remains a finite intensity 
at u = which transforms the gap to the pseudogap. The same behavior of the leading 
edge is observed in photoemission of underdoped cuprates which led to the idea of the 
photoemission pseudogap |l]-[| . 

In Fig. 8 we compare our calculated position of the spin-polaron peak and the ex- 
perimentally measured position of the photoemission leading edge as functions of 
momentum along the Fermi surface shown in the insets of both parts of the figure (some 
differences in the Fermi surfaces shown may be connected both with the experimental 
resolution which produces some uncertainty in the position of the Fermi surface espe- 
cially in the case of narrow bands and with the influence of terms not included in the 
t-J Hamiltonian). In part (b), the location on the Fermi surface is determined by the 
angle measured from the line (n, ir) — (tt, 0). At T = 14 K experimental curves for the 
83 K and 87 K samples in Fig. 8b correspond to superconducting gaps. As indicated in 
Ref . H , in the former, underdoped sample the shape of the curve and the magnitude of 
the gap remain practically unchanged when the temperature increases and somewhat 
exceeds T c , while the gap is closed in the latter, optimally doped sample. Taking into 
account these experimental facts we conclude from the comparison of Figs. 8a and 8b 
that the calculations reproduce satisfactorily the general shape and the magnitude of 
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Fig. 8. (a) The position of the spin-polaron peak along the Fermi surface, shown in the inset, 
for x = 0.121 (□) and x = 0.172 (O) at T = 116 K. The points M, X, and Y are defined 
in the caption of Fig. 6. (b) The position of the leading edge of the photoemission spectrum 
vs. momentum along the Fermi surface in the inset for samples with T c = 10 K (heavily 
underdoped, triangles), 83 K (underdoped, squares), and 87 K (optimally doped, circles). 
Measurements in Ref. ||, where this figure were taken from, were carried out for Bi2212 at 
T = 14 K. The points M and Y correspond to the X, Y and M points in the square Brillouin 
zone. 

the normal-state pseudogap in the underdoped case. Moreover, as seen from Fig. 8a, 
in agreement with experiment the pseudogap is closed in the case of optimal doping. 

In Ref. H the similarity of symmetries of the superconducting gap and the normal- 
state pseudogap and their smooth evolution into each other with temperature served 
as the basis for the supposition that the pseudogap is the normal-state precursor of 
the superconducting gap. In accord with this supposition the pseudogap arises above 
T c due to superconducting fluctuations. We did not include these fluctuations in our 
present calculations and thus the pseudogap shown in Fig. 8a is not connected with 
the fluctuations. In our calculations the pseudogap arises due to the specific dispersion 
of the spin-polaron band, a part of which is retained near the Fermi level at moderate 
doping and gives the most intensive maxima in the spectral function. The existence of 
this band is a consequence of strong electron correlations. We do not exclude the possi- 
bility that the superconducting fluctuations do contribute to the pseudogap, however, 
based on the satisfactory agreement between our estimate and the observed pseudo- 
gap magnitude we suppose that the main contribution is provided by strong electron 
correlations. Besides, with the superconducting fluctuation mechanism it is difficult to 
understand why the pseudogap appears only on one side of the Fermi level and why 
the fluctuations disappear abruptly in a small concentration range near optimal doping 
when the pseudogap is closed. For the 2D t-J model various calculations (see, e.g., 
Refs. p^ ,f£7|-|4"9|) give the d x 2_ y 2 symmetry for the superconducting gap in the case of 
moderate doping. This gap looks similar to the pseudogap in Fig. 8a, as the symme- 
tries of both of them are determined by the hole-magnon interaction and by short-range 
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antiferromagnetic ordering. This ensures smooth evolution of the pseudogap into the 
superconducting gap with lowering temperature below T c . 

As observed in Ref. [[|, in slightly underdoped samples the pseudogap is closed 
when temperature exceeds some T*. This characteristic temperature increases steeply 
with decreasing hole concentration. We believe that this behaviour is connected with 
crossing the boundary shown in Fig. 1. Outside of this boundary for x ^ 0.12 we found 
a hole energy spectrum which is similar to the spectrum of the usual metal and does 
not contain any pseudogap. 

Concluding this section let us mention two differences we see between our calculated 
and the experimental results. First, in our calculations the magnitude of the pseudogap 
increases slightly with decreasing hole concentration, while in the experiment || it 
remains unchanged, within experimental errors, for samples with markedly different x 
(see Fig. 8b). Second, in the normal-state photoemission spectra of underdoped Bi2212 
the linewidth increases dramatically on moving from (tt/2, 7r/2) to the vicinity of (0, it) 
@. Though the photoemission line shape is not well understood an analogous 
behavior of linewidths could be expected in the hole spectral function. However, in 
our results in Fig. 7 the change of the linewidth does not look so dramatic. This 
contradiction points either to the experimental pseudogap being somewhat larger than 
that reported in Ref. || [in this case the overdamped magnons wash away the (0, n) 
maximum, as they did with the (0, 0) maximum in Fig. 2a, curve 2] or to some decay 
process with low-frequency excitations not included in the t-J model. 



VII. MAGNETIC PROPERTIES 



The spin correlation function C\- 
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when 1 and m are on different sublattices. In the considered region of the T-x plane the 
decay of spin correlations (|19|) with the distance |1 — m| is nonexponential (see Fig. 9) 
which may be partly connected with finite-size effects. To estimate the spin correlation 
length we used the formula 

e 2 = Elil 2 ^ nl («>/(2Ee im («})- 



For T = 116 K the product $,^/x is nearly constant and approximately equal to the 
lattice spacing in the range 0.017 < x < 0.09, in agreement with experiment in 



La 2 _ :r Sr :E Cu04 [26]. For larger x at this temperature, £ becomes of the order of the 
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Fig. 9. The spin correlation function along the 
lattice axis [1 = (1,0)] for x = 0.027 (O) and 
x = 0.1 (□) at T = 0. 
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Fig. 10. The inverse correlation length in units 
of the inverse lattice spacing for x ~ 0.03 
(open circles) and x « 0.1 (filled circles). 
Values obtained from the neutron-scattering 
experiments in Lai.ggSrox-iCuC^ |5l| and in 
Lai.96Sro.o4Cu04 [§6| are shown by stars and 
diamonds, respectively. 

i i 

100 200 
T, K 

lattice spacing. As mentioned, the use of the spin-wave approximation (0) assumes 
that £ a. We hope, however, that our results obtained for the case £ ~ a (the 
region near the curve in Fig. 1 for x > 0.09) may give at least qualitatively a correct 
description of this region. When x ^ 0.06, our calculated £ is nearly independent of 
temperature in the considered range of T, as seen in Fig. 10 for the « 0.03 

(for fixed /i the hole concentration is somewhat changed with T; values of x given in 
the figure captions here and below are mean values for the considered temperature 
ranges). This behaviour agrees with experimental data in Lai. 96 Sr a04 CuO4 P6] , also 
shown in Fig. 10. For larger x values £ _1 are also weakly temperature-dependent at 
low T and increase more rapidly as the boundary in Fig. 1 is approached (see the data 
for x ~ 0.1 in Fig. 10). This behaviour is also in agreement with experimental results 
in Lai.86Sr .i4CuO4 by G. Aeppli et al. reported in Ref. |5l]] (stars in Fig. 10). Notice, 
however, that the growth of £ -1 predicted by the theory is somewhat more rapid than 
that observed in experiment. The saturation of £ with decreasing T, demonstrated by 
Fig. 10, is the distinctive property of the quantum disordered regime [H,2%| in which 
the system resides in the considered region of the T-x plane. 

In this region the hole contribution \ z h to the susceptibility is negligibly small in 
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Fig. 11. Imaginary x" an d real x' parts of the magnetic susceptibility as functions of wave 
vector along the symmetry lines in the Brillouin zone for u ~ 0.022 eV (a) and frequency 
for k = (vr,7r) (b). x ~ 0.11. Open and filled circles correspond to T = 58 K and 116 K, 
respectively. 



comparison with the spin contribution x z m - As seen m Fig- Hj * ne magnetic susceptibil- 
ity (|l~6D is strongly peaked around (7r,7r). This demonstrates strong antiferromagnetic 
fluctuations which persist even in the case of comparatively short correlation lengths 
of the order of the lattice spacing. As known, low-frequency incommensurate spin 
fluctuations are observed in the normal state of Sr-doped L^CuC^ fl52 |. This becomes 
apparent in a four-peaked structure of Imx(qu;), the peaks being displaced from the 
commensurate position to the points = (tt,tt±5), (ti ±5, it). To investigate whether 
this incommensurability is connected with the hole-magnon interaction in Cu02 planes 
a more sophisticated spin-wave approximation than that given in Eq. (|7]) is needed. We 
do not consider this point in the present paper. Notice that the frequency dependence 
of Imx(qu;) in Fig. lib is close to that observed in normal-state Lai.seSro.uCuC^ at 
wave vectors q« |51 |. Both the total peak intensity at these wave vectors, the position 



of the peak, and its temperature dependence are close to those shown in this figure. 

We used the obtained magnetic susceptibility for calculating the spin-lattice relax- 
ation times at the Cu and O sites 67 7i and 17 Ti and the Cu spin-echo decay time T 2 g 



by using the equations [14 
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F(q) = 2C 2 [1 + cosfe)]' 



where 63 7 n is the Cu nucleus gyromagnetic ratio, B = 3.82 ■ 10 7 eV, A± = 0.845, 



A\\ = —4B, and C = 0.915. The spin-lattice relaxation times in Eq. (p0|) correspond to 
the applied static magnetic field perpendicular to the CuC>2 plane. Results are shown in 
Fig. 12 together with the calculated x° — x(0, 0) and the respective experimental results 
obtained in YBa2Cu306.63 14|j53,|54 . From our results for different hole concentrations 



we selected for this figure those which appeared to be closest to the experimental data. 
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Fig. 12. The static spin susceptibility (a), 1/(T\T) at the Cu (b) and O (c) sites for H||c, and 
I/T2G (d) calculated for x 0.05 (open circles). The experimental data for YBa2Cu30g.63 
are shown by filled circles. x° was inferred from the Knight shift data H], 1/(T-\T) and 
I/T2G are from Refs. Q and |54j], respectively. 



21 



Those results are for x ~ 0.05 which is apparently somewhat smaller than the hole 
concentration in YBa 2 Cu 3 66 3. With increasing x our calculated x° an d ("TiT) -1 
increase, remaining of the same order of magnitude as those shown in Fig. 12, while 
T^q slightly decreases. These concentration dependences of x° an d T^q agree with 
experimental observations p!4| , p4| -|56|] , though for x° the theoretical dependence is 
stronger than that observed experimentally. The situation is apparently more diffi- 
cult for ( 63 TiT) _1 . In underdoped YE^CuaOe+y this quantity decreases with doping 



57J , while in HgBa 2 Cu04 +( 5, it increases p8 \. Apparently a more elaborate model is 



necessary for the description of this concentration dependence. Analyzing Fig. 12 we 
conclude that the t-J model is able to describe correctly the temperature dependences 
of the depicted quantities and to give their proper orders of magnitude in the quantum 
disordered regime. 

As seen in Fig. 12, (TxT) -1 and x° decrease with decreasing temperature. Such 
behaviour observed for T > T c in underdoped cuprates is considered as an indica- 



tion of the pseudogap in the spectrum of magnetic excitations @,|,|T^,|J]. Our results 
corroborate this point of view — analyzing Eq. and our numerical data we came to 
the conclusion that the mentioned temperature dependence is mainly connected with 
the occupation of low-frequency magnons, rather than the temperature variation of 
the magnon spectral intensities. Due to the pseudogap this occupation decreases with 
temperature. 



VIII. SUMMARY 

In this work we applied the modified spin-wave theory with the additional constraint 
of zero staggered magnetization to the two-dimensional t-J model in the paramagnetic 



state. In the Born approximation the constraint equation ( 12]) and the self-energy 
equations (|17D for the hole and magnon Green's functions form the self-consistent set 
which was solved numerically for finite hole concentrations x and temperatures T. The 
constraint can be fulfilled in the region of the T-x plane below the curve in Fig. 1. 

A number of unusual features of photoemission spectra in cuprate perovskites are 
satisfactorily reproduced by the obtained results. Among these features are the general 
shape of the electron energy spectrum and its evolution from the narrow spin-polaron 
band at small doping to the much wider band for moderate doping. Both by the 
energy position and by the extension in the Brillouin zone our calculated extended 
saddle point reproduces well the van Hove singularity of photoemission spectra. Also 
the obtained Fermi surface with hidden and two-dimensional parts is close to that 
observed experimentally. In the calculated Fermi surface the size of the hole pockets 
varies only slightly with the growth of the hole concentration. This growth is mainly 
connected with an increase of the quasiparticle weights of states occupied by holes. 
In our calculated hole spectrum the pseudogap has the magnitude, symmetry and the 
concentration dependence which are similar to those observed in photoemission. This 
pseudogap is not connected with superconducting fluctuations which were not included 
in the calculations. It arises due to the specific dispersion of the spin-polaron band a 
part of which is retained near the Fermi level at moderate doping and gives the most 
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intensive maxima in the spectral function. The persistence of this band is an indication 
of strong electron correlations retained at moderate doping. 

We calculated also a number of magnetic characteristics of the t-J model. We found 
that in the region shown in Fig. 1 the temperature variation of the spin correlation 
length is typical for the quantum disordered regime and that this dependence as well 
as the concentration dependence of the correlation length are close to those observed 
in cuprates. The magnon spectrum contains the pseudogap which manifests itself in 
the temperature dependence of the static spin susceptibility and spin-lattice relaxation 
rates at the Cu and O sites. Our calculated values of these quantities and the Cu spin- 
echo decay rate, their temperature and concentration dependences are in qualitative 
and in some cases in quantitative agreement with experiment. 

The considered phase with the above-discussed properties differs essentially from 
the conventional metal. As mentioned, this phase exists in the region bounded by the 
curve in Fig. 1. Outside of this region for x > 0.12 we found the energy spectrum of 
the usual metal. 
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